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ABSTRACT 

Context. Clusters are potentially powerful tools for cosmology provided their observed properties such as the Sunyaev- 
Zel'dovich (SZ) or X-ray signals can be translated into physical quantities like mass and temperature. Scaling relations 
are the appropriate mean to perform this translation. It is therefore, important to understand their evolution and their 
modifications with respect to the physics and to the underlying cosmology. 

Aims. In this spirit, we investigate the effect of dark energy on the X-ray and SZ scaling relations. The study is based 
on the first hydro-simulations of cluster formation for diferent models of dark energy. We present results for four dark 
energy models which differ from each other by their equations of state parameter, w. Namely, we use a cosmological 
constant model w = — 1 (as a reference), a perfect fluid with constant equation of state parameter w = —0.8 and one 
with w = —1.2 and a scalar field model (or quintessence) with varying w. 

Methods. We generate iV-body/hydrodynamic simulations that include radiative cooling with the public version of the 
Hydra code, modified to consider an arbitrary dark energy component. We produce cluster catalogues for the four 
models and derive the associated X-ray and SZ scaling relations. 

Results. We find that dark energy has little effect on scaling laws making it safe to use the ACDM scalings for conversion 
of observed quantities into temperature and masses. 
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1. Introduction 

In order to explain the current accelerati on of the universe 
(jRiess et al. 19 98; Pcrl mutter et al. 19991 ) in the context of 
the theory of General Relativity, it is general procedure to 
introduce a new form of gravitational component with neg- 
ative pressure - dark energy. Various candidates such as 
a cosmological constant or a quintessence field have been 
proposed. These models are characterised by their equation 
of state parameter and constrained by either using the ob- 
servation of background quantities or the growth of cosmic 
structures. An alternative way to explain the acceleration 
of the universe is to allow for modifications of gravity. Many 
classes of models exist, for instance, a light scalar field cou- 
pled to matter leads to models of extended quintessence 
and more generally to scalar-tensor type theories. Further 
possibilities were studied in the context of brancworld mod- 
els. Testing the Poisson equation on large scales may be a 
way to distinguish between all these alternative scenarios. 

Regardless of its nature, dark energy as a dominant 
component, plays a role in structure formation and thus 
modifies the number of formed structures. The evolution 
of linear perturbations in the presence of a scalar field like 
quintessence and the effects on the abundance of collapsed 
structures and its dependence with redshift were widely 



explored and suggested as a tool to constrain the nature 
of dark energy an d its e volution, e.g. jH aiman et al. 2001; 
I Weller et al. 2001t IWeinberg fc Kamionkowski 2003 
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iBattve fc Weller 2003t IWang et al. 2004t iMohr 20041) . The 
properties of collapsed halos (density contrast and virial 
radius) depend strongly on the shape of the potential, the 
initial conditions, the time evolution of the dark energy 
equation of sta te and on its ability t o collapse at the struc- 
ture scale, e.g. (|Nunes fc Mota 2 006). One can additionally 
investigate the effects of dark energy on the growth rate of 
structure and consequently study how dark energy affects 
the abundance of collapsed halos in dark en ergy models, 
e.g. (|Nunes et al. 20061 : iManera fc Mota 20061) . 

Cluster number counts can potentially discriminate be- 
tween dark energy models, in particular, the evolution of 
their equation of state. However, a large number of clus- 
ters with known redshifts from SZ, X-ray or optical surveys 
are needed, e.g. (|Bartelmann et alTl i005). Future Sunyaev- 
Zel'dovich (SZ) as well as X-rays observations will provide 
us with high quality data making galaxy clusters an effi- 
cient and powerful cosmological tool. Cluster physics, how- 
ever, is complex. The presence of substructures, the possible 
contamination by radio and IR sources, the still imperfect 
knowledge of the relation between the halo mass and the 
clusters observed properties induce degeneracies between 
cluster physics and cosmological models. Scaling relations 
are key quantities in observational cosmology as they re- 
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late the observables in X-rays and SZ to the cluster proper- 
ties, namely, masses and temperatures. The latter are then 
used in cluster number counts to constrain the cosmological 
models. 

Due to their importance in translating observations to 
physical quantities, but most of all in probing the cluster 
formation, scaling relations have raised considerable atten- 
tion. Simple models of formation of virialised systems such 
as clusters pr edict that the y exhibit self-similar b ehaviours 
dKaiser 1986h . see also (jAscasibar et al. 2006| ). In the 
self-similar model, gravitational infall drives shock heating 
of the intra-cluster medium (ICM) and establishes the gas 
properties such that they scale with the halo mass giving 
rise to the scaling relations. However, it is now clear that 
additional physics is required to provide a more complete 
picture of the cluster formation and evolution, and to 
explain the deviations between observations and the predic- 



tions based on self-similar scaling as 


shown mostly by X-ray 


observations (lEdge & Stewart 1991 


: Allen & Fabian 1998: 


Markevitch 1998 


Nevalaincn ct al. 2000: 


FinoEuenov et al. 20011 Ettori et al. 2004; iHenrv 2004: 


Arnaud et al. 20051: iRasia et al. 20051: iBaloeh et al 200a 


Maughan et al. 20061: Morandi et al. 20071). As the corn- 



plexity of the physical description increases due to the 
additional gas physics (galactic winds and/or quasar 
outflows, radiative cooling, preheating) the use of nu- 
merical simulations appears the best option to compare 
predictions and observations and examine the role of those 
new ingredien ts in explaining the departures form self sim- 
ilarity (e.g. (lEvrard et al. 1996t iBrvan fc Norman 1 998: 
Bialek et al. 20011: iThomas et al 20011; iBabul et al 2002; 
Voit et al. 2002t (Borgani et al. 200l iRowlev et al. 20041: 
Muanwong et al. 20061: iKav et al. 2007t )). The ICM can 
additionally be probed by the Sunyaev-Zel'dovich (SZ) 
effect, inverse Compton scattering of cosmic microwave 
background (CMB) photons off high energy electrons. 
The magnitude of the SZ effect is determined by the 
integrated gas pressure along the line-of-sight, called 
the Compton parameter, y. Unlike the X-ray surface 
brightness, the SZ effect is not subject to cosmological 
dimming and can be used out to high redshift which makes 
SZ scaling relations particularly interesting and attractive 
tests. Nevertheless, SZ measurements, although gaining in 
quality and quantity, are at the moment still not sufficient 
to be fully used and this explains why investigations of 
SZ scaling relations ar e still limited from the observa- 
tional point of view (ICoorav 19991: iBenson et al. 2004 
iMcCarthv et al. 20031 : iMorandi et al. 20071 ) as well as 
from 



the nume r ical point of view (|da Silva et al. 2004: 
IMolt et al. 20051 : IBonaldi et al. 20071 ). 

Scaling laws relate the observed properties in X-rays 
and SZ to the cluster properties, namely, masses and tem- 
peratures which are then used to construct mass functions 
and number counts utilised to constrain the cosmological 
models. Understanding the possible biases in the scaling 
relations is thus essential for the use of clusters as cosmo- 
logical probes. Previous studies, based on numerical simu- 
lations, have focused on the effects of additional gas physics 
in the scaling laws. In the present study, we explore the ef- 
fects of the cosmological model on the scaling properties of 
galaxy clusters and their evolutions. We perform hydrody- 
namic numerical simulation of cluster formation and evolu- 
tion assuming a simple radiative cooling model, rather than 
a more complete gas model, which allow us to single out 



only dark energy properties. We then investigate weather 
the X-ray and SZ scaling laws derived for different dark 
energy models depart from those obtained in the standard 
cosmological constant model (the ACDM model) taken as 
our reference. In the next section we briefly present the 
dark energy models used for the first numerical simulation 
of cluster formation with hydrodynamics in dark energy 
dominated universes. We describe the simulation code and 
the procedure used to construct the X-ray and SZ cluster 
catalogues. In Sect. 02 we present the X-ray and SZ scaling 
laws studied in the article together with the fitting proce- 
dure. Our results and conclusions are summarised in Sect. 
S] and \E\ respectively. 



2. Dark-energy simulations 

2.1. Simulation models 

Numerical N-body simulations including a dark energy 
component were performed by several groups to comple- 
ment the analytical computations of structure formation 
in presence of dark energy, and to the study the effects of 
dark energy at the structure level. All studies on galaxy 
clusters, were essentially dedicated to study dark mat- 
ter halo shapes and mass functions in different models 
of dark energy ([ Linder fc'jcnkins 2003; Lokas et al. 2003: 
iKlypin et al.2003t iDolag et al. 2004T~iKuhlen et al. 20041 ). 
The overall picture that has emerged from these stud- 
ies is that halo mass functi ons are well approx imated by 
the Jenkins mass function (|Jenkins et al. 20011 ) and halo 
core densities, or concentrations, are sensitive to the mean 
cosmological density at halo formation and therefore de- 
pend on the underlying dark-energy model. These find- 
ings have led to investigations of the effect of dark energy 
on the lensing properties of simulated cluster sized halos 
and their possible use t o constrain dark energy, see e.g. 
(|Meneghetti et al. 2005b1 : iMeneghetti et al. 2005at) . In the 
present study we perform the first hydrodynamics simu- 
lation of cluster formation in dark energy dominated uni- 
verse and we focus, for the first time, on the baryonic com- 
ponent of clusters and investigate the possible effects of 
dark energy on t heir gas properties. In an earlier study, 
(M aio et al. 2 006) produced numerical simulations of dark 
energy with baryonic gas to study the implications on cos- 
mic reionisation from first stars. 

In all the aforementioned studies, the scalar field associ- 
ated with dark energy is assumed not to have density fluc- 
tuations on scales of galaxy clusters or below. If dark energy 
influences th e perturbations on s mall scales as proposed fo r 
example by (lArbev et al. 20011). (IBean fc Magueiio 2 002) , 



(|Padmanabhan fc Choudhurv 2002t ) or (|Bagla et al. 20031 ). 
the collapse of structures itself will be affected. In our sim- 
ulations, we will also ignore any possibility for dark energy 
to cluster and influence the cluster formation. 

In addition to the cosmological constant model with 
a constant equation of state parameter w = pdc/Pdc = 
— 1, we simulate f our other models previously studied in 
(N unes et al. 20061 ). These phenomenological models span 
the range of values that the equation of state parame- 
ter can take for typical quintessence models and are com- 
pactible with current observational constraints. We take 
two models for which the dark energy is given by a 
perfect fluid with constant equation of state parameter: 
w = —0.8 and w = —1.2 (phantom dark energy), and 
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one model where dark energy results from a slowly evolv- 
ing scalar field in a potential with two exponential terms 
(2EXPl)( (jBarreiro et al. 200(f )) 



V{</>) = V (e Q ^ + e^) 



(1) 



where a = 6.2 and (3 = 0.1 which has equation of 
state parameter today wq = —0.95. One other model 
( 2 EXP 2) with var ying equation of state was studied in 
(|Nunes et al. 2 006). However, because its energy contribu- 
tion to the total energy of the universe quickly decays with 
redshift it has no significant departures from the Lambda 
model in the redshift range of cluster formation (e.g. z < 5). 
We thus choose, from this point onwards, not to consider 
this model it in this work. We further assume for all mod- 
els that energy density of dark energy, dark matter and 
baryons today are respectively fide = 0.7, £! m = 0.3, 
f2 B = 0.0486 and the Hubble parameter h = H /100 km s" 1 
Mpc- 1 = 0.7. cr 8 = 0.9. 

2.2. The simulation code 

For each of the cosmological models considered, we per- 
formed iV-body/hydrodynamic simulations of structure 
format ion using a modified v ersion of the public Hydr a 
code (jCouchman et al. 19951 : iPearce fe Couchman 19 97). 
which implements an adaptive particle-particle/particle- 
mesh (AP 3 M) al gorithm to calculate gravitational forces 
(jCouchman 199lh and smo othed particle hydrodynam- 
ics (SPH) (Monagha n 19921 ) to estimate hydrodynamic 
qua ntities. The SPH implemen tation follows that used 
by (|Thacker fc Couchman 2 000) and conserves both en- 
ergy and entropy. All simulations included a model 
for radiative cooling using the method described in 
([Thomas fc Couchman 199 3) and ba sed on the cooling ta- 
bles of ([Sutherland fc Dopita 1 993). The metalicity was 
assumed to be a global time varying quantity, Z = 
0.3(i/to)-Z©j where t/to is the age of the universe in units of 
the current time and Z@ is the solar metalicity. At a given 
time step, gas particles with temperatures below 1.2 x 10 4 K 
and overdensities (relative to the critical) larger than 10 4 , 
are converted into collisionless material and no longer par- 
ticipate in the gas dynamical processes. 

We modified the computation of the physical time, t, 
and scale factor, a, in Hydra to account for the effect of 
time- variable equations of state of dark energy, 



t = 



da' 



o a'H(a>) "° J E(a>) ' 
where E(a) = H(a)/H , or 



E 



( a) 2 = ^ +fidce - 3 />+^'»^ + 



isL 1 - f2 m - f2 



dc 



(3) 



for flat cosmologies. To speed up computations, this quan- 
tity is pre-tabulated for each of the dark energy models and 
read once at the beginning of the simulation run. E(a) is 
then interpolated and used in Eq. @. With these modifi- 
cations our version of the Hydra code can thus be used to 
simulate generic models of homogeneous dark energy. 

The initial density field of simulations was constructed, 
at redshift z = 49, using TV = 4, 096, 000 particles of bary- 
onic and dark matter, perturbed from a regular grid of fixed 



comoving size L — 100 h~ 1 Mpc. We use the Zel'dovich ap- 
proximation and the same set of random numbers to gen- 
erate the initial displacements. The amplitude of the mat- 
ter power spectrum was calculated assuming as = 0.9 in 
all models and the line ar growth factor comp uted for each 
model, as presented in ([Nunes fc Mota 20 06). We refer the 
reader to this article for a discussion on the normalisation. 
We also assume that the matter power spectrum transfer 
function is the same for all models and equals that of the 
cosmological co nstant model, which is computed using the 
BBKS formula (jBardeen et al. 19861) and the sh ape param- 
eter r given by the formula in ( Sugivama 19951 ). With this 
choice of parameters, the dark matter and baryon particle 
masses are 2.1 x 10 10 h~ 1 M Q and 2.6 x 10 9 h~ 1 M Q respec- 
tively. In physical units the gravitational softening was set 
fixed to 25ft, -1 kpc below z = 1 and above this redshift 
scaled as 50(1 + z) _1 h~ 1 kpc. 

Individual simulation runs took between 2592 to 2812 
time steps to evolve to z = 0. For each run we stored a total 
of 50 simulation snapshots (box outputs) at a list of selected 
redshifts (the same to all runs) ranging from z — 20 to z — 
0. Thirty of these outputs are inside the interval < z < 3, 
which is typically the range where galaxy clusters form. 



2.3. Catalogue construction 

From simulations, we constructed cluster catalogues 
using a modified version of the cluster extraction 
software devel oped at Sussex b y Thomas and col- 
laborators dThomas et al. 19981 : IPearce et al 20001 : 
iMuanwong et al. 2001 ). To summarize, the cluster 
identification process starts with the construction 
of a minimal-spanning tree of dark matter particles 
whose density exceeds the mean density of the box by 
Ab = 178 x f2~ - 55 (,z) (i.e., the density contrast predicted 
by the spherical collapse model of a virialised sphere 
relative to t he mean background densi ty in the Lambda 
cosmology ()Eke Navarro fc Frenk 1 998)). Although Ab 
may differ for different dark-energy models, this is not 
important at this step because cluster properties are 
computed at fixed overdensities as described below. The 
minimal-spanning tree is then split into clumps of particles 



— 1/3 

using a maximum linking length equal to 0.5 A b times 
the mean inter-particle separation. Finally we grow a 
sphere around the densest dark matter particle in each 
clump until the enclosed mass verifies 



(2) M A (< R A ) 



4tt 



R 3 A Ap crit (z), 



(4) 



where A is a fixed overdensity contrast, p C rit(z) = 
(3Hq /8ttG)E 2 (z) is the critical density and E(z) is given 
by Eq. ((3]). We have constructed master catalogues for all 
dark-energy simulation containing at least 500 particles of 
gas and dark matter, i.e. with an equivalent minimum mass 
of Mum « 1.18 x 10 13 h~ 1 M < 7 ) , at four fixed overdensities, 
A = 200, 500, 1000, 2500. Here we will report our findings 
for cluster scaling relations only at A = 200, the largest 
cluster overdensity radius usually considered in the litera- 
ture. For this catalogue we find 377, 393, 396, 374 clusters 
at z = in the cosmological constant, w = —0.8, 2EXP1 
and w = —1.2 simulation runs respectively. Although our 
choice of Ab may limit exact comparison of numbers, 
these abundances reproduce the behaviour predicted in 




Fig. 1. Cluster scaling relations Lx,200 — ?x,200 (left panel) and Y200 — -^200 (right panel) at redshift zero. Displayed 
quantities are computed within -R20O1 the radius where the mean cluster density is 200 times larger than the critical 
density. The embedded plots show the best fits with a power law to clusters represented in the main plots for the w = — 1 
(triangles), w = —0.8 (diamonds), 2EXP1 (squares) and w = —1.2 (circles) models. The shaded regions in the embedded 
plots give the typical scatter of the fits, i.e. the r.m.s dispersion around the best fit lines. 



(|Nunes et al. 200(1 ) in their analytical study. We note that 
the cluster d efinition used in this paper is different from 
that used in (jda Silva et al. 2004 iMuanwong et al. 20061) . 
Despite the similar simulation parameters, direct compar- 
ison with their findings is only possible for the Lambda 
model at redshift zero, where cluster definitions are identi- 
cal. 

For each model, the cluster catalogue provides us 
with estimated structural and observable quantities. More 
specifically for this study we compute: masses; intrinsic SZ 
luminosity, Y — Y$z X D\ where Ysz is the integrated 
SZ signal and £>a is the angular diameter distance; mass- 
weighted gas temperature T mw ; and bolometric X-ray tem- 
perature, Tx, and luminosity, Lx, excluding a cooling ra- 
dius of about 50 /t -1 kpc around th e cluster centre. We refer 
the reader to (jda Silva et al. 20041 ) for the definitions of the 
these quantities. 

3. Analysis of the scaling relations 

In this paper we investigate the following scaling relations 
between cluster properties: Tx — M, Y — M, Y — T mw , 
Lx — Tx, and Y — Lx- These can be expressed as: 

Tx = A TM (M/MoT™ (1 + zf™ E(z) 2 ^ , (5) 

Y = A YM (M/M a ) aYM (1 + zf™ E(z) 2/3 , (6) 

Y = A YT (T mw /r mw , ) QYT (1 + zf^ E(z)- 1 , (7) 
Lx = A LT (T x /Tx,o) aLT (1 + zf^ E{z) , (8) 

Y = A Yh (L x /ix,o) QYL (1 + z)^ E(z)' 9 / 4 , (9) 

where we have chosen the normalisation scales Mq = 
lO^/i^M®, T x ,o = T mw , = 1 keV, L x ,o = 10 43 



erg/s/h 2 . The powers of the E(z) function give the pre- 
dicted evolut ion, extrapolated from the self-similar model, 
(jKaiser 19 86). of the scalings in each case. The quantities, 
A, a, and /?, give: the scalings normalisation at z = 0; the 
power on the independent variable; and the departures from 
the expected redshift evolution. 

To investigate these cluster scaling relations in 
our simulations w e use the method described in 
(jda Silva et al. 20041 ). According to Eqs. (j5ll9"j) . the general 
form of how a given cluster property y relates to a property 
x can be written as, 

yf(z)=y (z)(x/x ) a , (10) 
where 

y (z) = A(l + zf, (11) 

and f(z) is some fixed power of the cosmological factor 
E(z). This is a power-law function whose parameters A, 
a and (3 can be obtained by fitting our cluster catalogue 
distributions at each redshift with a straight line in the 
(log(y f(z)),\ogx) plane. To be more specific, the fitting 
procedure is carried out in three steps. Firstly, we fit the 
cluster distributions with a straight-line in logarithmic scale 
at all redshifts. If the logarithmic slope a remains approxi- 
mately constant (i.e. shows no systematic variations) within 
the redshift range of interest, we then take a at z = as 
the best fit value. In the second step, we repeat the fit us- 
ing this value of a to determine the scaling normalisation 
factors yo{z)- This avoids unwanted correlations between a 
and yo(z). Finally, in the last step we use Eq. (fTTj) to obtain 
the parameters A and [3. 

In the fitting process we consider only clusters with 
M\im > 5 x 1O 13 /i _1 M0 for scalings with mass, Lu m > 
6.6 x 10 42 /?,~ 2 erg s _1 for scalings with luminosity, and 
Imw.iim > 1 keV, Tx,iim > 1-1 keV for scalings with mass- 
weighted and emission- weighted temperatures, respectively. 
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This ensures that our cluster samples are complete at all 
redshifts and that, for each model, equal number of clus- 
ters are used for all scalings at redshift zero. With these 
selection criteria our cluster samples contain 60 clusters for 
the cosmological constant model, and a similar number of 
clusters for the other models, at z=0. We note that above 
z = 1.5 the number of clusters with Mu m > 5 x 10 13 ft- _1 M Q 
in our sample decreases typically below 10, hence, we do 
not fit the scaling relations above this redshift value. As we 
will discuss in the next section, all the scaling relations ex- 
plored in the present study, are well fitted by power- laws of 
the form Eq. (fTO)) . in the redshift range < z < 1.5, except 
for the Lx — 2x and Y — Lx whose normalization depen- 
dences with redshift, \og(yo(zj), are well approximated by 
a straight line only in a narrower redshift range. 

4. Results 

4.1. Scaling relations at z = 

We start by discussing the cluster scaling laws at redshift 
zero. These were determined, for all scalings and models 
under investigation in this work, at an overdensity radius 
i?200- In Table Q] we present the power law best fit for the 
logarithmic slopes a obtained for each of the cases. As can 
be inferred from the table, all models provide very similar 
results for each scaling, with differences between models be- 
ing comparable or in most cases within the statistical errors 
of the fit. This indicates that the local cluster scaling re- 
lations are quite insensitive to the underlying dark energy 
model driving the present-day evolution of the universe. As 
stated in Sec. 12.31 our cluster definition is the same as in 
(|da Silva et al. 2004 iMuanwong et al. 20061) (only) at red- 
shift zero, for the Lambda model. We verified that indeed 
theirs and our results are in excellent agreement at this 
redshift. 

To illustrate the robustness of the scalings with respect 
to the dark energy models investigated in this paper, we 
present in Figure Q] two characteristic X-ray and SZ galaxy 
cluster scaling relations: the Lx — Tx (left panel) and the 
Y — M (right panel) scalings. In each case, the main plot 
shows the cluster distributions for all models whereas the 
embedded panels show the power law best fits obtained. For 
both scaling relations, the cluster distributions and best fit 
lines clearly overlap. Also represented by a shaded area in 
the embedded panels is the r.m.s. dispersion of the fit for 
Lambda model: 



(12) 



where y' = yf (see Eq. (fl"U|) ') and y[ are individual data 
points. This dispersion is of the same size of the fit disper- 
sions obtained in the other models and it is clearly wider 
than the best fit line separations of the various dark en- 
ergy models. As expected, the scatter in the (core excised) 
Lx — Tx is larger than in the Y — M relation due to the 
higher sensitivity of the former scaling to the gas physics 
in the inner regions of clusters. 

4.2. Evolution of the scaling relations 

In Figs. [H [31 and [5] we present our findings for the vari- 
ation with redshift of a (top panels) and yo(z) (middle pan- 
els) for the T x - M, Y - M, Y - T mw , L x - T x and Y - L x 
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Fig. 2. Slope (upper panel), normalisation \og(yo(z)) (mid- 
dle panel), and normalisation best fit lines (bottom panel) 
of the Tx— M relation as a function of redshift (\og(yo(z)) is 
defined in Eq. ITTj). Blue colour and triangles stand for the 
cosmological constant model, cyan and diamonds are for 
the w = —0.8 model, green and squares are for the 2EXP1 
model, and yellow and circles are for the w = —1.2 model. 
The shaded area in the bottom panel gives the dispersion of 
the normalisation fit for the cosmological constant model. 



relations, respectively. Here, triangles, diamonds, squares 
and circles represent the w = —1, w = —0.8, 2EXP1, and 
w = —1.2 models, respectively. The coloured bands give the 
the best fit errors obtained at each redshift for these quan- 
tities. The lines in the bottom panels of these figures are 
linear best fits to the evolution of \og(yo(z)) with log(l + z). 
The shaded area in these panels gives the rms dispersion of 
the \og(yo(z)) fit for the cosmological constant model (simi- 
lar values of dispersion are found for the other models). The 
resulting best fit parameters for, A, (3 and a are presented 
in Table [TJ for all scaling investigated in this paper. 

As can be inferred from the figures, the fit to a power 
law, in a redshift range < z < 1.5, was possible for all 
the scaling relations explored in the present study except 
the relations involving Lx (namely Lx — Tx and Y — Lx)- 
For the latter two, we found no significant departures of the 
slopes a from the cosmological constant model but we could 
not fit the evolution of the scaling relations by a power law 
within the whole redshift range < z < 1.5 (see Figs. [5]- [5]). 
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lines of the Y — M relation as a function of redshift. The 
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For these scalings we therefore restrict the analysis to the 
a smaller range, namely < z < 0.75, ie (log(l + z) < 0.25) 
where the linear fit is valid. 

We compared the average value of the slope a of each 
scaling relation, for all the cosmological models, over the 
redshift range < z < 1.5 and found no significant depar- 
tures from one model to the other, confirming the behavior 
observed at z = 0. This is also clear from the upper pan- 
els of Figs. [3J [31 dj E] and O where there's a high degree of 
overlap between the a's obtained for the four cosmological 
models. For scalings with mass (T x — M and Y — M) the 
degree of overlap at high redshift is somewhat less striking 
due to larger variations caused by the rapid decrease of the 
number clusters with M\[ m > 5 x 10 13 . These "oscillations" 
show however no systematic dependences with redshift and 
are about the same "mean values" for all models. 

In order to study the evolution of the scaling relations, 
we factored out the redshift dependence expected from the 
self-similar evolution which is parametrised by a power law 
with exponent (3, see Eq. (jlip . As a general remark, we 
find that, excluding the T x — M relation, all scalings in- 
vestigated here show positive evolutions relative to the ex- 
pected self-similar evolution (i.e. for a given x in Eq (flO)) 
the property yf is higher at higher redshift). In the case of 



T x — M relation the opposite behaviour is observed. This 
is expected because the inclusion in all models of radia- 
tive cooling (which is a non-gravitational physical process) 
causes cluster scaling laws to deviate from self-similar evo - 
lution, see eg (|da Silva et al. 2004 iMuanwong et al. 2 006) 
for studies in cosmological constant model simulations. 

It is clear from the table and the lower panels of the 
figures that the value of (3 is slightly more model dependent 
than the parameters a and A. This is especially the case 
for the Lx — Tx, and Y — Lx relations for which there is a 
mild difference between models. For these scalings the w=- 
1.2 and w=-0.8 model generally show the largest deviations 
from the cosmological constant model, whereas the 2EXP1 
show the smallest. However, those differences are of the 
same order as the intrinsic errors and dispersions. We can 
thus safely consider that there are no significant departures 
from the cosmological constant model. 

5. Conclusions 

The abundance of clusters, their redshift distribution, as 
well as their clustering, are governed by the geometry of 
the universe and the power spectrum of the initial density 
perturbations. Gas physics related to cluster structure and 
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evolution also enters through mapping of the cluster observ- 
able (SZ flux or X-ray luminosity) relative to the total mass 
of the cluster. As a result, galaxy cluster counts can be used 
as probes of cosmological and cluster properties. However, 
there are several requirements needed to achieve precise 
cosmological constraints: (i) advances in understanding the 
formation and evolution of cluster size halos; (ii) a good 
understanding of the selection function; (iii) robust obser- 
vational proxies for the cluster mass. 

The first condition relates to the conduction of large 
simulations. This is greatly achieved in the standard cosmo- 
logical model with a cosmological constant. In the context 
of dark energy dominated universe, N-body simulations are 
becoming available for models different from the simple cos- 
mological constant with constant or varying equation of 
state parameter. These simulations now provide us with a 
good understanding of the halo properties. They show that 
the halo mass function is well approximated by the Jenkins 
mass function. Simulations also indicate that the clusters 
halos are more concentrated in dark energy models, since 
structure grow earlier, than in lambda models and that 
concentrations are higher in models with varying than con- 
stant equation state. The second condition, understanding 
the cluster selection function, translates in understanding 



the limiting mass and the completeness of the surveys from 
realistic mock cluster catalogues. The mass (or tempera- 
ture) selection function is directly linked with the cluster 
observed quantities through cluster scaling relations which 
is our third requirement (the need for a good proxy for the 
cluster mass). 

In this work we have thus for the first time explored 
the scaling laws for both SZ and X-rays observations using 
hydrodynamic simulations of galaxy clusters in four dark 
energy models with constant or varying equation of states 
spanning a large class of models. We have studied the scal- 
ing properties at z = and their evolution with redshift. 
We have found that dark energy induces no modifications 
on the scaling laws at z = and presents very little dif- 
ferences from the cosmological constant model at higher 
redshifts. 

While detailed simulations incorporating viable dark en- 
ergy models remains a program in progress, it is reassuring 
that all models considered in this work predict similar scal- 
ing properties to the lambda model. The modeling of the 
cluster gas component appears to be nearly independent 
of the dark energy model. Therefore, using the "standard" 
lambda model scaling relations for converting observable 
to masses and temperature in future surveys should not in- 
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troduce any additional bias in the cosmological constraints 
derived from cluster counts. In this work, we have consid- 
ered that dark energy does not cluster with dark matter. It 
would be interesting, however, to evaluate how our conclu- 
sions stand for numerical simulations in a scenario where 
dark energy is inhomogeneous and collapses along with dark 
matter during the formation of structure. This will be pur- 
sued in a forthcoming study. 
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Table 1. Best fit values of the parameters a, log A and (3 as well as their respective la errors. These values are valid 
within the rcdshift range < z < 1.5. For the Lx — Tx and Y — Lx relations, the linear fit and the associated parameters 
are only valid in the range < z < 0.75 above which the linear fit is not a good approximation. 
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Ptm 


0.620 ± 0.029 
0.260 ± 0.005 
-0.228 ± 0.020 


0.604 ± 0.031 
0.267 ± 0.004 
-0.249 ± 0.017 


0.581 ±0.033 
0.271 ±0.005 
-0.264 ±0.023 


0.602 ± 0.029 
0.258 ±0.005 
-0.230 ± 0.023 


Y - 


-M 


ttYM 
log A Y M 

Pym 


1.732 ±0.025 
-5.910 ±0.004 
0.128 ±0.016 


1.730 ± 0.025 
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0.116 ±0.016 
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log A Y t 
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2.922 ±0.100 
-6.522 ± 0.008 
0.454 ± 0.036 
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-6.518 ±0.005 
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2.838 ± 0.072 
-6.499 ± 0.007 
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Lx 




QLT 

log Alt 
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2.738 ± 0.086 
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2.589 ±0.011 
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1.063 ±0.028 
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1.064 ± 0.026 
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0.890 ± 0.028 
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1.084 ± 0.037 
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0.497 ± 0.037 



